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I. INTRODUCTION 


At least two experimental evidences demand for physics beyond the standard model 
(SM). The first one comes from the well stablished neutrino oscillation experiments [1-4] 
which imply that all the three known neutrinos (v e , z/ M , u T ) are quantum superpositions 
of three massive states zy (* = 1, 2, 3). The second evidence is firmly stablished from 
several observations and studies of gravitational effects on different scales, which points out 
that most of the Universe’s mass consists of non-baryonic dark matter (DM). Specifically, 
the Planck collaboration has determined that the DM relic abundance is given Ddm^ 2 = 
0.1193 ±0.0014 [5], 

In order to explain both of these evidences, it is now clear that the SM has to be extended. 
In the neutrino case, usually new fermionic fields, Ur s, are introduced to generate Dirac mass 
terms for neutrinos. The Ur fields are, in general, singlet under the SM gauge groups, and 
thus these can also have Majorana mass terms. Moreover, in order to explain the smallness 
of active neutrino masses, the Ur s usually get large masses via the well known see-saw 
mechanism [6, 7]. On the other hand, the existence of DM in the Universe requires at least a 
new massive particle since SM does not provide any viable DM candidate. The most studied 
and well motivated candidates for DM are the weakly interacting massive particles (WIMPs). 
In general, WIMPs are neutral, stable and are present in a plethora of extensions of the SM 
[8-11]. Nowadays, there are several astroparticle experiments actively pursuing detection of 
WIMP DM candidates in direct and indirect ways. The direct detection experiments [12— 
14] have set upper bounds on WIMP-nucleon elastic scattering, whereas the indirect ones 
[15-19] presented upper limits on the thermal average of the same scattering cross section 

In this work we study a scenario that simultaneously offers an explanation for the previ¬ 
ously mentioned open questions of the SM. In special, we present a local B — C extension 
of the SM in which there are three Ur fermionic fields and some extra scalars, 4>, 0*. The 
B — C quantum numbers of the extra new fermionic fields come from exotic solutions of 
anomaly constraints. These solutions were found for the first time in Ref. [20]. Here, we 
propose a simplified version of the model in Ref. [20, 21] where an almost natural Z 2 sym¬ 
metry stabilizes ur 3 and thus allows it to be a DM candidate. Another appealing feature 
of this model is that it implements a see-saw mechanism at low energy because neutrino 
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masses are proportional to V£/V<p (where and are VEVs of the <L, 0*, respectively) 
and Vp has been set in MeV energy scale. The discrete symmetry also simplifies the task 
of setting the Yukawa couplings in the neutrino mass Lagrangian in order to agree with the 
neutrino oscillation parameters. There have been extensive studies on these two matters, in 
special in gauge extensions of the SM (a few of which are contained in the Refs. [22-24]). 
This is so because U(l) gauge factors are contained in grand unification theories [25, 26], 
supersymmetric models [27] and left-right models [28, 29]. 

The paper is organized as follows. We start by discussing the model in Sec. II. In that 
section we present its held content and general Lagrangian. We also show the almost natural 
Z ,2 symmetry in the model which stabilizes the DM candidate n^R- In Sec. Ill we study 
the scalar sector in detail. We obtain analytical formulas for both the mass eigenstates and 
the eigenvalues when it is possible. We also include a discussion about the consequences of 
the presence of the Majoron J which is due to the breaking of an accidental global U(l)j 
symmetry. Specifically, we show that it escapes the current bounds on energy loss in stars 
[30, 31], effective number of neutrinos N e g [5], and the invisible decay widths of Higgs [32-35] 
and Z\ gauge boson. In Sec. IV we analytically find the parameters of the neutrino mass 
matrices in order to satisfy the data from the neutrino oscillation [36] and other constraints 
such as lepton flavor violation (LFV) [37, 38]; the sum of the SM neutrinos masses [5]; and 
effective Majorana mass m ee from double beta decay experiments [37, 38]. In Sec. V we 
carry out a study of the relic dark matter abundance and the direct detection prospects. A 
general discussion follows in Sec. VI where we present our conclusions. Finally, in Appendix 
A we show the general minimization conditions coming from the scalar potential. 

II. THE MODEL 

We consider an extension of the SM based on the gauge symmetry 
SU(2) L <g>U(l) y /<g>U(l) B _£ where B and C are the usual baryonic and leptonic num¬ 
bers, respectively, and Y' is a new charge different from the hypercharge Y of the SM. The 
values of Y' are chosen to obtain the hypercharge Y through the relation Y = \Y' + B — C ], 
after the first spontaneous symmetry breaking. The fields of this model with their respective 
charges are shown in Table I. Actually, this model is a simplified variation of the one 
introduced in Refs. [20, 21]. Specifically, here we have removed one of the extra doublets 
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of scalars considered there. As we will show below, this allows an almost automatic Z 2 
symmetry that stabilizes the DM candidate, n^. The remaining scalar fields are enough to 
give mass to the neutrinos at tree level. It is also important to note that there is an exotic 
charge assignment for the B — C charges where (B — C-) nR1) n R2 — ~4 and (B — C) n , l3 = 5 
different from the usual one where ( B — C) nm = 1 with i — 1, 2, 3. 
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Table I: Quantum number assignment for the fields in the model. I 3 , Y' and B — C are the quantum 
numbers under the symmetry groups SU(2)^, U(l) y /, U(1 ),b-£, respectively. 


With the field content in Table I, we can write respectively the most general renormaliz- 
able Yukawa Lagrangian and scalar potential respecting the gauge invariance as follows 

--2V = Yf l) L Li e Ri H + Y<f>Q Li d Rj H + Y^Q Li u Rj H + V im L Li n Rm $ 

+^M mn (7i Rm ) c n Rn (p 1 + -M 3 3 {n R3 ) c n R3 (p 2 + ^M m 3 {n Rm ) c n R3 (f) 3 + H.c., (1) 

and 

V B -c = + X H \H^H\ 2 — |$ f $| 2 — /i 2 a |0„| 2 + A q |0*0 q | 2 

+k h * \H \ 2 |$| 2 + + n Ha \H\ 2 |0 a | 2 + |$| 2 |0 O | 2 

+«a^(0Q0a) (000/?) + [«1230102(03) 2 _ + &123X (0X01) (0203) 

+ K 3x(0X03) + H.c.] , (2) 

where i , j = 1, 2, 3 are lcpton/quark family numbers; m, n — 1, 2; H — (t 2 is the Pauli 

matrix), and a,/3 — 1,2,3,X with a ^ j3 in the /t a/ 3(0Q0a:) (0^0,3) terms. Also, we have 
omitted summation symbols over repeated indices. 
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Before we go further, two important remarks are in order. Firstly, from Eqs. (1) and 
(2) we see that apart from the \M. m z{ n Rm) c nRs^ + H.c. terms, the Lagrangian is invariant 
under a Z 2 symmetry acting in a non-trivial way on the n r% held, i.e. Z 2 (n^ 3 ) = —Ur 3 (the 
rest of holds being invariant under this symmetry). We will consider the case of this Z 2 
symmetry throughout this work. Hence, the Ur 3 fermionic held will be the DM candidate. 
Secondly, from Eq. (1) we see that quarks and charged leptons obtain masses just from the 
H vacuum expectation value, (H°) = Vr. Therefore, the H interactions with quarks and 
charged leptons are diagonalized by the same matrices as the corresponding mass matrices. 
In this case the neutral interactions are diagonal in havor and there is no flavor-changing 
neutral current in the quark and charged lepton sector. This feature remains after the 
symmetry basis is changed to mass basis [39, 40]. However, lepton havor violation (LFV) 
processes coming from the terms proportional to T>i m and T> i:i can occur at one loop. We 
will discuss these processes in more detail in Sec. IV. 

III. SCALAR SECTOR 

In the general case this model has a rich scalar spectrum and its vacuum structure can 
take several configurations. However, we are going to make some simplifying and reasonable 
assumptions that allow us, in most cases, to obtain analytical formulas in both the neutrino 
and the dark matter sectors. We will discuss systematically our assumptions throughout 
this paper. 

Firstly, as result of the absence of one of the extra doublets and of writing only the 
renormalizable terms in the scalar potential, the model here considered has a Majoron, J, 
in its scalar spectrum. This is a general conclusion and does not depend on any particular 
choice of the set of parameters. Once the neutral scalars develop non-vanishing vacuum 
expectation values, VEVs, and using the usual shifting tp () = (V^ + Rey? + i Im </?) for the 

scalar fields (the superscript “ 0 ” means we are taking the neutral part of the held), we End 
that J can be written as 


J = ^- 9V2V H V£e 2 Im H° - 9 V2V%V£e Im 4>° 



Im 4> 2 


+3V2V 2 ( V 2 + V 2 e 2 ) Im </> 3 + 9y/2V* (V 2 + V 2 e 2 ) lmcp x 


(3) 
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where Nj = V 2 ^J(AV 2 + (3V 2 + 4I0 2 ) e 2 ) (58 V* + (3V 2 + 58V|) e 2 ) and e = 10/10. We 
also have defined the VEVs as (00 = V^ 7 with z = //,$, 0i, 0 2 , 03, 0x and set V0 = I0 2 = 
10, 3 = I0 X = 10 for simplicity. The parameter e is chosen <C 1 as we will show below. 

The presence of J in the physical spectrum is due to an extra symmetry in the scalar 
potential in Eq. (2). In other words, the scalar potential actually has a larger global 
SU(2)i<g)U(l)y/<g)U(l)B_£<g)U(l)j symmetry. The last symmetry, U(l)j, acts on the scalar 


fields H , <3>, 0i, 02 , 03, 0x with charges — 23 ’ 23 , ff, respectively. We have nor¬ 

malized the charges in order to set the 0i charge equal to 1 . Also, note that the global 
symmetry U(l)j is independent on the U(l) y ' and U(1 )b_£ symmetries. This is neces¬ 
sary to consider it as an actual extra symmetry. Furthermore, U(l)j can be extended 
to the total Lagrangian acting on the fermions Ql, u R , d R , L l , e R , n Rm , n 3R with charges 
0, —1|, ||, —1|, —0 — A, respectively. Therefore, J is a true Majoron with mass equal 

to zero at all orders in perturbation theory. Gravitational effects can break this symmetry, 
and thus give mass to the Majoron [41-45]. Studies taking into account these effects on 
B — C symmetry constrain the energy scale of its breakdown to be < 10 TeV [44], However, 
we are not going to consider this case. 

The major challenge to models with a Majoron comes from the energy loss in stars through 
the process 7 + e~ —» e~ + J. This process is used to put limits on the eeJ coupling, g Se j, 
and it is found that it must be g Se j < 10 “ 10 for the Sun, and g Se j < 10 “ 12 for the red-giant 
stars [30, 31]. In our case, g Se j = —- 777 -^ e 2 where Y0 ' 1 and m e are the 

electron Yukawa coupling to the H scalar and electron mass, respectively. Since e = 10/10, 
10 m = >/ Vh + Vi and V H ^ V su (the R eH° is the only field giving mass to the top quark 
at tree level), we have that t«l. Thus, expanding g SeJ in series of e, it is straightforward 
to see that g Se j — 2 v 7 ^y >2 + 0(e 4 ). Choosing 10 = 1 TeV and V R — Vsm = 246 GeV we can 

notice that e < 3.8 x 10 -4 in order to satisfy the limit coming from red-giant stars analysis. 
It is straightforward to show that the smallness of e is technically natural, since doing e —>■ 0 
increases the symmetry of the total Lagrangian. 

The charged sector can also be found analytically. Besides the charged Nambu-Goldstone 
eaten by the W ± gauge boson, the model has one charged scalar, C ± . It can be written as 

(J± — _I _ H± _j_ T/„ (Tli _ KfiMpfl _j_ K 'h<S> V R 


KHi>xV£ 


, .0 n = (10e E ± + V H < h ± ), with squared mass given by m 2 c ± = + 

\/W+R> e v 1 e 


+ \k' h §V£ e 2 . Note that when e —* 0 we have that in general rric± —> 00 . However, 
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when this happens, the minimization conditions in Appendix A require that Kh^x oc e. 
Thus, mc± remains finite. 

In order to find the rest of the mass eigenvalues and eigenstates of the scalar potential 
(the CP—even, CP—odd scalars), in general, we numerically proceed choosing the set of 
the parameters to satisfy simultaneously the minimization conditions given in Eqs. (Al- 
A6), the positivity of the squared masses, and the lower boundedness of the scalar potential. 
All these constraints are always checked numerically. Furthermore, we restrict ourselves to 
a relevant set of parameters that allows us to study the dark matter properties in some 
interesting cases. Our initial assumptions are: (i) For the sake of simplicity: = V^ 2 = 

V<t > 3 — V<t> x = ( we have already used this in Eqs. (3) and in the C ± charged scalar), 

Km = = K H1 = KH3 = KhX = = ^<S> 2 = = K$X = K 12 = ^ 1.3 = KlX = K 23 = 

hi 2 x = « 3 x = 0 and H 123 X = n' 3X — ^123 (h) In order to have the heaviest CP—even scalars 
with similar masses, we choose: Ai = A 2 = A 3 = Ax = Aand (iii) Due to the stability 
of the minima, we obtain: Hh<s>x = V^,e (see Eq. (A2)) and /P H = A hVh + ~ e 2 ; 

A 4 ! = -+ ^'J , ^ 2e2 > A 4 ? = (Pl23 + p 2 = H 2 2 H + ( K 123 + ^<t> )^0> A 4 3 = (3^123 + ^0)^, 

and fi 2 x = (K 123 + A <p)V£ — -|^e 2 . The rest of parameters will be set when required. 

In general, the squared mass matrices of the CP—odd scalars (M^, p _ odd ) and the 
CP—even scalars (M^ p _ even ) can be written in powers of e up to e 2 , i.e. Mf = + 

eMjA + e 2 M|j with i = CP — odd, CP — even. In spite of the smallness of e and the as¬ 
sumptions made above, it is a hard task to obtain exact analytical expressions for the mass 
eigenvalues and mass eigenstates of these matrices. These can be found perturbatively in 
powers of e, though expressions are usually very long and no more clarifying. In this section 
we just provide the leading-order expression of the scalar masses because these yield a good 
picture of their exact behavior. 

In the CP—odd sector the model has three scalars, C, / 2 , I 3 , besides the Majoron J in 
Eq. (3) and the two Nambu-Goldstone eaten by the Z\ (it is assumed that Z\ is the gauge 
boson with mass equal to the Z boson in the SM) and Z 2 boson. Their masses are given by 
rrir t = Sj/V ^ j q ’ , mi 2 = a/5 — V7y/—n 123 V 0 , rrij 3 = ^5 + ^/7 s/—k\taVo- From the previous 
expressions we see that k \23 < 0 in order to have all masses belonging to reals. It is also 
straightforward to see that I\ = Im$° + (9(e). Additionally, we find that J 2 and J 3 are, at e 
order, a linear combination of the Inin’s with i = 1,2,3, X. The CP—even sector is more 
complicated even in the leading order. In this sector the model has six different eigenstates, 
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Ri s, with masses given by: m ^ 1 = \/2X hVr, tur 2 = , rURs = V 2 ^<t> ~ 3 - 58 1^123 1 V^, 

m RA = a/2A^ + 1.15 [K 123 TV 0 , m ^ 5 = y/2 (A^ + |/ci23|)V^, rriR 6 = %]2X^ + 2.42 | ^ 123 1 ^ / 0 - -Ri 
(which is Re 77° + O(e)) is the scalar that plays the role of the Higgs scalar boson in this 
model, since it couples at tree level to all fermions, giving mass to them when it gains a 
VEV, Vff. Thus, we set its mass equal to 125 GeV. We find that Xr — 0.13 — 0.14 gives 
the correct value for the Higgs mass. R -2 is Re$° + 0(e). The rest of fields are, in general, 
combinations of the Reel’s with i = 1, 2, 3, X. Note that the remaining four CP— even 
scalars have masses proportional to V/, and there is not a criterium to determine precisely 
their masses. However, we have to choose the parameters in the scalar potential such that 
all rriR masses are larger than the 2 b boson mass (m% ~ - r~^ 2 ir^ = "I'n ) due to the 2 b 
invisible decay width. In other words, if some of mR i were < M Zl then the Z\ boson could 
decay through the process Z\ —> Ri + J —> J + J + J, which would contribute to the Z\ 
boson decay width as half of the decay Z\ —» vv [46]. According to the experimental data 
there is no room for such an extra contribution [36]. 

All expressions above for masses and eigenstates are very useful to have a general view of 
the scalar spectrum. However, it is necessary to work with more precision when calculations 
of the DM sector are involved. Thus, from here on, we always work numerically to diagonalize 
the squared-mass matrices for both the CP—odd and the CP—even scalars. 

Finally, a further comment regarding the J presence is necessary. Since the Majoron 
J is massless, it contributes to the density of radiation in the Universe which is usually 
parameterized by the effective neutrino number N e g. This parameter specifies the energy 
density of relativistic species in terms of the neutrino temperature. Planck together with 
WMAP9 polarization data, high-/ experiments and the BAO data (Planck + WP + highL 
+ BAO) gives N e s = 3.301q|i [5]. In the case that the Majoron J goes out of equilibrium 
when the only massive particles left are electrons and positrons it makes a contribution to 
A r ef f equal to 4/7 which is in agreement with the current data. In the case when J decouples 
in higher temperatures a lower contribution is expected. For a best treatment see Ref. [47] . 











IV. NEUTRINO MASSES 


The mass Lagrangian for neutrinos, which conies from Eq. (1) when the neutral scalars 
gain VEVs, can be written in matrix form as: 


_ = - 


r -[ 


0 

M D 


1 - 

O 

VL (■ n c ) L 








[Ml 



n R 


H.c., 


( 4 ) 


where vl = [u ]y n , n R = [ri\ n 2 n 3 ]^. The Majorana and Dirac mass matrices ( M m 
and Mn-, respectively) are written as 


"" = ^ 


^ M u M 12 0 ^ 


■Ml2 Ai 22 0 

0 0 M-33 


nr V * 

MD = T 2 


( 


\ 


'D ii V 12 0 

T*21 T^22 0 

T )31 ^32 0 


\ 


( 5 ) 


For e 1, i.e. <C V^, the mass matrix in Eq. (4) can be diagonalized using the regular 
see-saw mechanism. The masses of the heavy neutrinos, iVj with i — 1,2, DM, are related to 
the energy scale of the VEVs of the singlets and are given by the eigenvalues of Mm- Mn 2 = 
(M n + M 22 ) =F v / 4Af ? 2 + (^ii-^ 2 2 ) 2 ] V (2y/2) , M Na = M^/y/2 = M DM - For 
simplicity, we set Ml 12 = 0 and Ain = Ai 22 . Doing so, we have M Nl2 = M N = AinV^/y/2 
and M N3 = Mdm- We work with Mdm , Mjy and as input parameters. Thus, Ain 
and Af 33 are expressed in terms of Mdm , Mat and as Ain = V^Mn/V^ and AI 33 = 
\/ 2 Mdm/V^, respectively. 

As it is well known, the masses of the light neutrinos, zy with i = 1 , 2 , 3, are given by the 
eigenvalues of the matrix M v = M D M^M^. From Eq. (5), it can be seen that det M v = 0. 
It implies that at least one of the light neutrino masses is zero. The minimal requirement 
for the parameters in AIm and Ain is that these have to provide the light neutrino masses 
and mixing angles consistent with the oscillation neutrino constraints. There are other 
constraints on neutrino masses such as < 0.23 eV coming from Planck collaboration 

[5] that we are going to consider below. 

Now, we proceed analytically making the ansatz that M u is diagonalized by the tri- 
bimaximal-Cabbibo (TBC) matrix, I/tbc [ 6 ], i.e. U^ bc M u Utbc — M u = diag(mi, m 2 , m 3 ). 
For the sake of simplicity, we parametrize mi = x — y, m 2 = 2 x + y, m 3 = 2 1 / + x — y, and 
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work in a basis where the charged lepton mass matrix is diagonal. C/tbc can be written as 


/ 


^TBC — 


1-4 

3 \ 1 4 


'7h ( 1 + a ) 


JL h-4 

V3 l 1 4 


-4 \ 


A 

75 


14 4 

x/3 V 1 2) ^2 l 1 4 


1_ 

^ 7 e ( 1_A ) + 72 ( 1 _ t)J 


4 

4 

A 2 

4 


+ 0( A 3 


( 6 ) 


where we have chosen 5 = 0 (see [ 6 ]). 1 /tbc leads to the mixings: sin 0 12 = sin 6*13 = 
and sin 023 = We choose A = 0.218174 to be consistent with the experimental limits 
[36] for neutrinos. It is remarkable that this A value is consistent with the relationship 
sin 9c ~ A where is the Wolfcnstein parameter. I 11 order to set the parameters in Mm and 
Mo, we match and M D M^f as it should be if our ansatz is supposed 

to work. Note that we have used UtbcU^bc = 1 + 0( A 4 ) and Utbc^tbc = 1 + 0( A 4 ). In 
addition, we have one more degree of freedom to choose because the neutrino mass hierarchy 
is yet unknown. The neutrino mass hierarchy can be either normal (mi < m 2 < m 3 ) or 
inverted (m 3 < m\ < m 2 ). We separately consider them. 

In the case of normal hierarchy we choose m j = 0 since det M u = 0. Doing so, m 2 = 3 y 
and m 3 = 2v. Hence, we find 


UtbcM u U tbc — 

16 


( y (—4 + A 2 ) + A 2 z/ 2 (-4 + A 2 ) (y (-2 + A) - 2Xu) 
* (—2 + A) (4 y + (2 + A) u) 


2 (-4 +A 2 ) ( y (2 + A) — 2\v) ^ 
(—4 + A 2 ) (Ay + (—4 + A 2 ) 12 ) , 

(2 + A ) 2 (4 y + (—2 + A ) 2 v) j 


(7) 


where means that the assigned matrix element is equal to its transpose element. The 
matrix M D M^ is written as 




= K 


7A?i + 7A? 9 7 Aii7A 2 i + D 12 X? 22 


\ 


* 

* 


D 2 i + D, 


2 

22 


A>2lT>3i + V 22 P32 


7^31 + ^32 


/ 


( 8 ) 


y2 V 2 e 2 

where we have defined the dimensional constant K = t - = y(y= Matching Eq. (7) 

to Eq. ( 8 ) we have a system of six independent equations. We were not able to solve 
analytically that system for the six general variables X>n, 7A 12 , D 2 i, V 22 , 7A 31 , D 32 . However, 
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if we set T> 12 = 0 4 , we can solve it analytically for the remaining five variables and obtain 
the following four solutions: 


T>n — Si x 



(9) 


T> 2 i — Si x 


4 yfK' 

(A 2 — 4) ((A — 2)y — 2Xu) 


( 10 ) 


2 yjKr 


TX >22 — S 2 x 


\J (A 4 — 16A + 16) 2 uy 2 


( 11 ) 


V 3 i — si x 


Ay/Kry 

(A 2 — 4) ((A + 2 )y — 2A v) 
2 y/Kr 


( 12 ) 


^32 — S 2 X 


(A 4 + 16A + 16) yj (A 4 - 16A + 16) uy 2 


(13) 


Ay/Kry (A 4 - 16A + 16) 


where Si = { — 1, —1, +1, +1}, s 2 = {—1, +1, —1, +1} and r = 16X 2 u + (A 2 — A) 2 y. 

Now that we have the solutions for the T> 13 in terms of y, v and K , let’s find y and v from 
Am 2 un = m 2 — m 2 = 9 y 2 and |Am 2 tm | = m| — m 2 = Au 2 , with A?n 2 mi = 7.53 x 10~ 5 eV 2 and 
|Am 2 tm | = 2.52 x 10 -3 eV 2 [36]. We solve these equations obtaining y ~ 2.89252 x 10” 3 eV 
and v ps 2.50998 x 10~ 2 eV. Neither y nor v may be negative because m 2 and m 3 are positive. 
We then find the masses values: m\ = 0, m 2 ~ 8.677556 x 10“ 3 eV, m 3 ~ 5.01996 x 10“ 2 eV 
which shows a normal mass hierarchy. 

In order to determine completely the T) tJ values, we still have to find K. From Eq. (8) 
and assuming 0(T>ij ) ~ 1, we have that K < 10~ 10 GeV sets the neutrino masses in sub-eV 
mass scale. However, the value of K can not be taken arbitrarily small due to one-loop 
induced processes violating lepton flavor (LFV). Specifically, we consider LFV processes 
such as li —> lj + 7 , where i — fj,, t and j = e, /i, respectively. This model has one-loop 
contributions to these kinds of processes since charged leptons couple to charged scalars 
and right-handed heavy neutrinos. The branching ratio is estimated as Br (/* —>• lj + 7 ) = 



upper bounds for Br (/i — > e + 7 ) and Br (r —>• y + 7 ) are < 5.7 x 10 13 and < 4.4 x 10 8 


[36], respectively. The dependence of the Br (/j — » l 3 + 7 ) on K value arises through f M1 
and J'ei which depend on T>, 3 values (see Eqs. (9-13)). Also, note that Br (^ — > lj + 7 ) 


1 We have numerically solved the equations finding always 2?i 2 = 0, what justifies our choice. 
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weakly depends on Mjy k because F 2 (x) does not drastically depend on its argument. Thus, 
we can find a lower bound for K imposing the experimental upper bounds Br (// — y e + 7 ) 
and Br (r —>■ p + 7 ). We find that for 200 GeV < M Nl , 2 < 1000 GeV, the K value has to 
be > 2.72 x 10~ 9 GeV. We have also used mc± = 416.986 GeV, which is the value used 
in the DM analysis and is its correct limit at (9(e). Now, we can finally find the values for 
Dij. For instance, for the first solution displayed in Eqs. (9-13), we have Du ~ —0.0754, 
D 2 i ~ —0.1492, V 2 2 ~ —0.125, D :il ps —0.042 and V 32 ~ —0.1947. ft is important to 
say that the lower limit on K imposes a constraint on e > 1.97 x 10 -6 (we have fixed 
Mjv 12 = 1000/a/ 2 GeV and V^=l TeV). This constraint on e and the one coming from the 
safety of the Majoron J (see Section III) imply that 1.97 x 10 -6 <f<3.8x 10” 4 . 

Experiments on Ozy/3/3 constrain the effective Majorana mass \m ee \ = 
| c i 3 ('micf 2 e tSl + m^s^e* 52 ) V m 3 e 2l ^ CP s\ 3 \ and the strongest one is up to now 
\m ee \ < 0.27 — 0.65 eV, 90% C. L. [37, 38]. We are not interested in CP violation 
nor phases in the leptonic mixing matrix, therefore we end up finding m ee = 0.00376994 eV. 
Also, data coming from Planck collaboration [5] constrain the sum of the light neutrinos to 
be lower than < 0.23 eV. Clearly, it is satisfied since m\ + m 2 + m 3 0{ 10- 2 ) eV. 

For the inverted mass hierarchy, the procedure is very similar to the one shown above, 
thus, we present only the main results. In this case, we have m 3 = 0, as opposed to m\ = 0. 
It yields the parametrization rrii = x — y, m 2 = 2 x + y and m 3 = 0. For the D^ values, we 
have D 12 = 0 and 


D\\ — si x 


(A 2 - 4) V® 


D22 — D 3 2 — S2 x 


4 VK 

V( x ~ y)( x + 2 v ) 


D21 — Si X 
H31 = Si X 


Xx — 2 y 
2 VKx ’ 
Xx + 2 y 


(14) 


(15) 


V2A% ’ M 1 2 VK~x ’ 

where si, S 2 and K are defined as above. The values for x and y are found from Am 2 un = 

ml — ml = Sy(2x + y) >0, and |A??r, 2 tm | = m| — m\ = (x — y) 2 . If x > y, we have 
3y(2x + y) = 7.53 x 10“ 5 eV 2 and x — y — (2.52 x 10' 3 ) 1 / 2 eV. Solving for x and y , 
we have x ~ 5.04478 x 10~ 2 eV and y ~ 2.48162 x 10~ 4 eV (there is another solution: 
x ~ 1.64850 x 10~ 2 eV and y ~ —3.37146 x 10~ 2 eV, but we choose the first one to work with). 
We can then find the mass values: mi ph 5.01996 x 10 -2 eV, m 2 ~ 5.09441 x 10^ 2 eV, m 3 = 0 
which show an inverted mass hierarchy, as said in the beginning. The LFV processes require 
that K >1.77 x 10 9 GeV. Using these values for K, x, y and A, we obtain, for the 
first solution displayed in Eqs. (14-15), D n ~ 0.0523, V 12 ~ 0, V 21 ~ —0.0056, V 22 ~ 
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—0.0379, T >31 ~ —0.0061 and H 32 ~ 0.0379. Similarly to the normal hierarchy, we have that 
1.58 x 1CT 6 < e < 3.8 x 1CT 4 in order to satisfy the experimental bounds for Br (/1 — > e + 7 ) 
and Br (r —>■ g + 7 ). 

For the m ee limit, we hnd m ee = 0.0492258 eV, which is below the latest experimental 
limit [37, 38]. Regarding the Planck limit [5], it is satisfied because mi + m 2 + m 3 ~ 0.101 
eV. 

V. DARK MATTER 

As previously mentioned, this model has an almost automatic Z 2 symmetry acting on Ur 3 , 
i.e. Z 2 (iir 3 ) = —iiR 3 . We have imposed it to be exact in the total Lagrangian by removing 
just one term. Thus, Ur 3 is stable and it can, in principle, be a DM candidate. From here 
on, we consider IVdm (which is equal to % 3 , the difference being that N-£, M is a mass basis 
field and the former a symmetry basis one) as a DM candidate and verify whether it satisfies 
the current experimental data. These data come essentially from investigations of Planck 
collaboration [5] which constrains the DM relic density to be Ddm ^ 2 = 0.1193 ± 0.0014; and 
from direct detection (DD) limits of LUX [12], XENON100 [13] and SuperCDMS [14], which 
constrain the cross section, for scattering off nucleon, to be smaller than 7.6 x 1CU 10 pb for 
WIMP mass of 33 GeV. We will consider these constraints below. 

A. Relic Abundance 

In order to hnd the present DM relic density, Ddm^ 2 , coming from the IVdm Majorana 
fermion, we must solve the Boltzmann differential equation. This standard procedure is well 
described in Refs. [49, 50]. Here, we are not going to enter in its details because we have 
used the packages Feynrules [51], Calchep [52] and MicrOMEGAs [53]. The first two being 
auxiliary to the third that calculates Ddm ^ 2 for a given model which contains WIMPs. 

In Fig. (1), we show the processes which mainly contribute to the DM annihilation cross 
section, and so lead to the present relic density. All of them depend on the parameters 
in the Lagrangians given in Eqs. (1), (2) and on the kinetic terms involving the covariant 
derivatives. We have already fixed most of those parameters in Secs. Ill and IV. However, 
g, gY'-, gB-L, A h , AA$, hri, M 23 and Mdm remain still free. The first three parameters 
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g, gy-, 9b-l are the gauge coupling constants of the SU(2 )l, U (l) y / and U(1 )b~c groups, 
respectively. Roughly speaking, these couplings and the VEVs together determine the masses 
of the gauge bosons. The VEVs have already been set in the previous sections. In addition, 
g can be set equal to 0.652 due to the W ± mass, gy and gB-L mainly determine the 
masses of the Z 2 gauge boson and its mixing with Z\ in the neutral current. From precision 
electroweak studies [54-56], its mixing, given by tan/3, has to be < 10 -3 (see Ref. [22] for 
an analytical expression of tan/3). Furthermore, M Z 2 /gB-L > 6 TeV [57, 58]. We find that 
working with g Y ' = 0.506 and gs-L = 0.505, we obtain tan /3 ~ 2 x 10 -4 and M Z2 — 4.7TeV, 
as well as the known SM gauge bosons masses. 

Now, the A h parameter is chosen to be 0.13 < A h < 0.14 because it is the main responsible 
for the Higgs mass, Mm ggs = 125 GeV, when kh 2 < 0.1. In principle, the A$ and A</, 
parameters can take a wide range of values. Here, we have set A$ = 0.5 and A</, = 0.8, 
and thus we have the non-SM scalar masses larger than the SM particle masses. The kh 2 , 
K 123 and M dm parameters have been scanned in a broad region of values. Specifically, 
we have iterated the MicrOMEGAs package for the DM mass within the range 10 GeV < 
Mdm < 1000 GeV, taking into account different values of k,h 2 and K 123 , and leaving the 
remaining parameters constant. In general, we have worked with k,h 2 = 0.1, 10~ 2 , 10 -4 
and —0.56 x A^ < K 123 ^ 0 (0.56 x 0.8 = 0.448). The last choice because we must assure 
that all the scalar masses are real (we obtain a slightly more constraining condition on k 12 3 
if we impose that all CP—even scalar must have masses larger than the Higgs boson, i.e. 

(' m Hig g s — ^A^V^) ~ —0.442 < K 123 V 0). Also, it is important to note that fiq 2 3 
controls the scalar trilinear vertices between scalars. 

Regarding the Bh 2 parameter, we find that, in our scenario, it largely governs the in¬ 
visible Higgs width T|^g gs to non-SM particles. It is because k # 2 induces mixing between 
Re H° and Re </> 2 and thus it mostly determines the coupling Higgs- JJ, Chjj, since J has a 
component in Im 0 2 . This Chjj coupling induces a tree-level contribution to the T^ ggs given 
by C^jj/327miHigg S - Under the assumption that the production and decays of the Higgs 
are correctly described by the SM aside perhaps from decay into new unobserved particles, 
the branching ratio for the Higgs decay into new invisible particles, Br^ gs , is known to be 
< 10 % — 15% [32-35]. As kh 2 < 0.2 we find that the Br]j ggs remains under this value for 
—0.442 < K 123 ^ 0. We have been conservative choosing k H2 <0.1 for all results. 

Taking into account all aforementioned considerations on the parameters, we plot, in 
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Figure 1: Feynman diagrams which represent the main annihilation processes that contribute to 
the calculation of floM h 2 . We have defined: x = ^dm; Si = Ri, /*; A fl = Zi^, ; and A ° 

Fig. 2, Q DM h 2 versus M dm for k # 2 = 10” 1 , 10” 4 , with n 12 3 = —0.4 (the figure on the 
left) and k 123 = —0.1 (the figure on the right), respectively. The region in gray is ruled 
out because Ddm ^ 2 is overabundant. The dot-dashed line is the Ddm ^ 2 = 0.1193 value 
reported by Planck. In general, we find that depending on the Mdm, various annihilation 
channels are important and clearly some resonances are visible. Resonances are, in general, 
found in m mediator /2. Thus, for convenience, we give here the scalar masses for both figures 
in Fig. (2). For the case with ft 123 = —0.4 (both values of k , R2 ) we approximately have 
m Ri ~ 125.0, 417.0, 411.3, 1435.8, 1549.2, 1603.0 GeV, M h ~ 417.0, 970.4, 1748.8 GeV, 
and Mc± — 417.0 GeV. On the other hand, for the case with K 123 = —0.1 (both values 
of «ff 2 ), we approximately have rn Rt ~ 125.0, 417.0, 1114.6, 1309.7, 1341.6, 1357.3 GeV, 
Mj. ~ 417.0, 485.2, 874.4 GeV, M$± ~ 417.0 GeV. In all cases we have the Majoron J. 

In order to better comprehend the annihilation processes and their contributions con¬ 
tained in the curves in Fig. 2, we plot Fig. 3 which shows the relative contributions to 
Odm h 2 of the main DM annihilation channels. Let’s consider some relevant regions. For 
Mdm less than 80 GeV we have in general two resonances. The first one is due to the in¬ 
terchange of the Z\ gauge boson in the s— channel. It is located at Mdm = M^/2 ~ 45.6 
GeV and remains there even when Kui is set to 10 -4 . It is so because it depends on the 
TVdm — jV D m — Z\ coupling via neutral currents. Since this coupling arises from the covariant 
derivatives, it is independent on the Kh 2 value. In contrast, the second resonance, which 
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Relic density vs. DM mass (/C 12 3 =-0.4) 


Relic density vs. DM mass (*123 =-0.1) 




Mq m (GeV) 


Mq m (GeV) 


Figure 2: Figures displaying the dependence of Ddm ^ 2 on Mdm, for K 123 = —0.4 (the figure on 
the left) and K 123 = —0.1 (the figure on the right). Each figure shows two cases corresponding 
to kh 2 = 10 1 ,10 -4 . The dot-dashed line is the Planck QdmA 2 . The gray region means the DM 
candidate would be excluded, and the white one means that the DM candidate can still make some 
part of the DM content of the Universe. 
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Figure 3: (Color online) Figures displaying the relative contributions (%) of the main annihilation 
channels to the DM relic abundance. The cases with kh 2 = 10 _1 , 10 ~ 4 and K 123 = —0.1, —0.4 are 
shown. 
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arises by the s —channel interchange of the Higgs boson (located in mmggs/2 ~ 62.5 GeV), 
disappears when k H2 = 10 -4 . This is understood realizing that iV DM couples to the Higgs 
boson via the term 1 ^dm (nR 3 ) c nji 3 (j) 2 and since Higgs component in 0 2 depends on k, H2 , 
it is clear that the smaller kh 2 the smaller AA> m — ^dm - Higgs coupling. In this region of 
masses we also notice that Adm annihilation processes into quark-antiquark pair (in special 
into bb quarks for kh 2 = 10 _1 ) are the dominant for both figures. These occur via Higgs 
mediation. Annihilation processes into neutrinos via the Z\ interchange are also important 
(~ 25%). This is also true for both values of k H2 and for both k 123 = —0.4 and —0.1. 

As M dm increases from 80 GeV to 120 GeV, and as long as k i23 = —0.4 and k H2 = 
10 _1 , the annihilation into gauge bosons {W ± /Zi) are dominant (in particular into W + W ~) 
with some considerable (~ 20%) contribution of annihilation into quark-antiquark pair. In 
contrast, for n\ 23 = —0.4 and Kh 2 = 10~ 4 , the Adm annihilation processes into quark- 
antiquark pairs continue being the most important. Moreover, Adm annihilation processes 
into JJ start to be considerable (~ 15%). Similar conclusions are true for the case of 
^123 = —0.1 and kh 2 — 10 _1 . However, for this case, annihilations into gauge bosons 
(W ± /Z\) have a little lower contribution when compared to the case K\ 23 = —0.4. For the 
case of k 323 = —0.1 and k,h 2 = 10 -4 , the annihilations into antiquarks-quarks are the most 
contributing. 

In the region 120 GeV < M dm < 180 GeV and with k 123 = —0.4 and k,h 2 — 10 _1 , 
roughly speaking, three Nn m annihilation processes are similarly predominant. These are 
annihilations into W + W~/Z\Z\, R\R\ and JJ. Recall that R\ is the Higgs-like scalar. For 
this region of mass and with fiq 23 = —0.1 and k,h 2 = 10 _1 analogous conclusions can be 
reached. This is not the case for Kh 2 = 10 -4 (with fiq 23 = —0.4) in the same Mdm region, in 
which case, IVdm annihilations into JJ are almost completely dominant with an additional 
contribution (~ 12%) from the annihilations into quark-antiquark pairs. For k 123 = —0.1 
and kh 2 — 10“ 4 , however, we find that annihilation into quark-antiquark pairs are still 
dominant. 

As Mdm is around mn 3 /2 ~ 205 GeV, we see a resonance, in Fig. 2 on the left, due 
i? 3 s—channel interchange for both kh 2 values. Here, predominant annihilation process is 
A^dm-^dm —>• JJ with more than 50% of contribution. It is also important to note that for a 
M D m in this region we have the Udm ^ 2 value reported by Planck. This resonance does not 
occur in the k 123 = —0.1 cases because rrifi 3 ~ 1114.61 GeV. However, in these cases, we 
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have one resonance at mj 2 /2 « 242.6 GeV, with IVdm-Adm —» R\J as dominant process for 
Kh 2 = 10 _1 , and with N DM N DM —> qq as dominant process, for k H2 = 10~ 4 . 

In the region 220 GeV < M dm < 500 GeV and with k 12 3 = —0.4 and kh 2 = 10' 1 , we 
can say that two annihilations, Vdm-Vdm —>■ J J and -/Vdm-Adm —f P3J, strongly control 
I4 dm^ 2 - Except when Mdm ~ m / 2 /2 ~ 485.1 GeV where -ZVdm-Adm —>■ P 3 J annihilation 
completely governs Udm^ 2 - As 500 GeV < Mdm < 700 GeV, annihilations into JJ, R\I \, 
P 3 P 3 , R 3 J and JR are predominant and their contributions depend on the proximity to the 
three different resonances. Finally, when 700 GeV < Mdm < 1000 GeV, annihilations into 
JJ, R 3 I 11 R 3 J, N 1 N 1 and N 2 N 2 are the most contributing processes to determine m^ 2 - 
Similar behavior is found for the case K 123 = —0.4 and kh 2 = 10~ 4 . It is so because in these 
regions of masses the annihilation processes mostly depend on the trilinear vertices between 
scalars. 

As K 123 = — 0.1 and kh 2 = 10 _1 , the scalar spectrum changes and thus the location of 
the resonances change as well. As it was commented, the resonance at rrif i;i /2 « 205 GeV 
does not exist anymore. Instead, we have a resonance at mj 2 /2 « 242.6 GeV. In the region 
220 GeV < Mdm < 470 GeV, the most important difference, in contrast with the case 
AC 123 = —0.4, is that we do not have regions with Udm ^ 2 < 0.119 (a little tiny region can be 
seen in the mj 2 /2 « 242.6 GeV). Another difference is that annihilation into W + W~/Z\Z\ 
remains to be important in this region (~ 15% — 35%). In addition, annihilation into R,\I\ 
contributes > 35% in most of this mass region. Other annihilations, such as R\ J, I\Ii, J J 
and R\R\, also contribute but are subdominant. For 470 GeV < M D m < 700 GeV, Odm^ 2 is 
completely determined by annihilation into JlJl, hJ, JJ. As 700 GeV < Mdm < 1000 GeV, 
annihilation processes into N\N\ and N 2 N 2 share importance with RIi, R 3 I 1 , R 4 I 1 and R 5 I 1 
to determine Odm^ 2 - For Mdm > 470 GeV, we have Odm ^ 2 < 0.119. 

Finally, when ^23 = —0.1 and K 112 — 10” 4 , we have some relevant differences. What 
is clearest is that for M D m < 470 GeV we just have the Z x resonance, which depends only 
on the VEVs and the g' s. This is because of the smallness of K 123 and kh 2 (specially 
K-H 2 = 10~ 4 ), which makes the couplings with these CP—odd scalar mediators be tiny. Some 
features are worth mentioning, though. Up to Mdm — 350 GeV, annihilation into quarks are 
predominant. After that, until M D m — 700 GeV, the final products RIi, I\ J, JJ (summing 
~ 35 — 100%) enter as the major contributors to the relic density and the quarks enter as 
subdominant processes fading out at M D m — 450 GeV. Next, up to M D m — 900 GeV, the 
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Figure 4: (Color online) 2-D figure displaying the behavior of Ddm/i 2 as one continuously varies 
both K 123 in the range [—0.4, 0], with kh 2 = 0.1, and Mdm hr the range [10, 1000] GeV. The cyan 
points represent correct relic density Ddm^ 2 , within experimental errors; the green ones indicate 
Ddm ^ 2 below the Planck result; and the gray ones mean Ddm ^ 2 above Planck constraint and thus 
ruled out. 

main annihilation products are N\N\ and N 2 N 2 (~ 30 — 40% each), with R 3 I 1 taking place 
at the end of this interval. Finally, for 850 GeV < Mdm < 1000 GeV, the main contributions 
come from R 3 I 1 , R 4 I 1 and R-Ri, summing more than 70% of the DM annihilation energy. 

Now, in order to grasp the behavior of the relic density when one varies K 123 , we show 
a two-dimensional figure, Fig. 4, which was obtained with MicrOMEGAs, from a 10 5 points 
iteration. We see from it that as one varies k 12 3 , the regions for correct relic density (cyan 
points), change place, getting to the minimal value of Mdm ~ 200 GeV for K 123 = —0.4; and 
also for some points which increases in K 123 as Mdm decreases, having at K 123 ~ —0.05 its 
last point. We can also notice green regions (together with cyan lines) that extend from left 
to right as M D m increases, and the reason behind it are the resonances of J 2 , 1 3 (decrease as 
K 123 increases) and R 3 (increases). Therefore, one can conclude that the correct relic density, 
before Mdm ~ 500 — 600 GeV, may only be reached through resonances of the lightest singlet 
particles of our spectrum. 
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Scaled cross section vs. DM mass (A'- )2 3=-0.4) 
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Figure 5: Figures displaying the curves representing the SI cross section per nucleon, crsi-nucieom 
asa function of Mdm, for the IVdm elastic scattering off nucleon. Cases for K 123 = —0.4; —0.1 and 
kh 2 = 10 _ 1 ,10 ~ 2 are shown. In these figures, we also display the SI upper limits coming from 
LUX (dashed), XENONIOO (dot-dashed) and SuperCDMS (dotted). All <rsi_ n ucieon curves in the 
gray region are ruled out by the LUX upper limit. 

B. Direct Detection 


Other important constraints on DM candidates come from the current experiments [12-14] 
which aim to directly detect WIMP dark matter by measuring the kinetic energy transferred 
to a nucleus after it scatters off a DM particle. All of these experiments have imposed 
limits on the WIMP scattering cross section off the nuclei. In general, the WIMP-nucleus 
interactions can be either spin-independent (SI) or spin-dependent (SD). Currently, the most 
constraining limits come from the Large Underground Xenon (LUX) experiment [12] which 
has set bounds on the SI WIMP-nucleon elastic scattering with a minimum upper limit on 
the cross section of 7.6 x 10 -10 pb at a WIMP mass of 33 GeV/c 2 . 

We have verified that, for IVdm considered here, the dominant interactions are SI. Thus, 
we calculate (using the MicrOMEGAs package) the SI elastic scattering cross section per 
nucleon, Osi-nucieon, and the results are shown in Fig. 5. Actually, we scale the Osi-nucieon 
cross sections with the calculated relic density relative to that measured by the Planck in 
order to properly compare the predicted cross sections with those given by direct detection 
experiments, which present their results assuming the observed density. The experimental 
limits on SI cross sections are also shown in Fig. 5. We have not shown results for SD cross 
sections because we found that those are, in general, several orders under the SI current 
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limits, see Refs. [59, 60] which state a minimum upper bound of ~ 5 x 10 _3 pb at a WIMP 
mass of 24 GeV/c 2 . 

From Fig. 5, it can be seen that the smaller n H2 , the smaller Osi-nucieon- As Kr 2 — 10~ 2 
the a si-nucieon is below the LUX upper bound for all values of Mdm- For i%h 2 = 1CT 1 and 
Mdm < 500 GeV, crsi-nucieon is below the LUX limit only around the resonances. In contrast, 
for Mdm > 500 GeV, the LUX limits are satisfied for all cases shown in Fig. 5. This 
implies that Osi-nucieon mainly depends on n H2 . This fact is easily understood by realizing 
that, in our case, the relevant interactions for direct detection are mostly mediated via 
Higgs in the t-channel . Thus, these interactions depend on the mixings between the Higgs 
scalar (Ri) and rest of R{ scalars. These mixings strongly depend on the k,h 2 value, as it 
was already discussed. In addition, we can see from Fig. 5 that although asi-nucieon (actually 
ex si-nucieon x Ddm^ 2 /0.1193 ) does not depend directly on other scalars, there is clearly indirect 
dependence on them because these scalars affect the annihilation cross section and thus the 
relic abundance. 

Finally, from Figs. 2 and 5 we can conclude that, provided nH 2 1CU 2 , the constraints 
coming from the Hdm^ 2 determine if a set of parameters leads to a viable dark matter 
candidate or not. 

VI. CONCLUSIONS 

In this paper, we have discussed a scenario where neutrino masses and dark matter are 
possible. In particular, the model presented here is a gauge extension of the SM based 
on SU(2 ) l <g) U(1)y/ ® U(l) B _ jC symmetry group. Besides the SM fields, we have added 
one doublet scalar (<F), four singlet scalars (0j), and three right-handed neutrinos ( n Ri ). 
These three last fields have different quantum numbers under the gauge groups and more 
importantly, they couple to different scalars. This allows a rich texture in the neutrino mass 
matrices. In addition, because of the exclusion of one term from the Lagrangian, we have 
that a Z 2 symmetry acting only in urs appears. This opens the possibility that rira be a 
DM candidate. 

The model contains a very rich scalar sector, and in special, we find that it contains a 
Majoron, J, which has its origin in a global accidental symmetry, U(l)j. We show that 
this symmetry is exact and it acts on the total Lagrangian. Despite the fact that global 
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symmetries can be broken by gravity quantum effects, we have not considered that possibility 
in this paper. Thus, the Majoron remains massless and in principle, it poses some issues 
to the safety of the model. Therefore, we study the consequences of the presence of J in 
the physical spectrum. Specifically, we consider four major challenges: (i) energy loss in 
stars by the process 7 + e~ —> e~ + J [30, 31]; (ii) relativistic degrees of freedom in the 
Universe, parametrized by N e g] (iii) Z invisible decay width; and (iv) Higgs invisible decay 
width. Setting — 1 TeV and V R — Vsm = 246 GeV, the first of these constraints leads to 
e = V^/Vfj) < 3.8 x 10 -4 (where we have set all VEVs of scalar singlets equal). The second 
constraint really does not lead to any restriction on the parameters of the model since the 
J contribution to the density of radiation in the Universe is in agreement with the Planck 
limits [5] for different decoupling temperatures. The issue of the Z invisible width is also 
overcome by imposing some constraints on the scalar potential parameters as shown in Secs. 
Ill and V. In addition, the Higgs invisible decay width is maintained at Br]j" gs < 10 — 15% 
[32-35] provided that k,h 2 < 0.2. Finally, since e < 1, we have shown an overview of the 
scalar spectra by expanding the scalar squared-mass matrices in powers of e and doing other 
simplifying considerations. 

Since e has to be very small, we can use the well-known see-saw approximation to analyt¬ 
ically solve the T>ij and A4 t j parameters. These are found by imposing some experimental 
constraints coming from the neutrino physics. In special, the mixing angles and the differ¬ 
ences of the squared neutrino masses. We manage to solve the and A iij parameters by 

making the ansatz that M v matrix is diagonalized by the tri-bimaximal-Cabbibo matrix. It 
is important to note that the existence of the Z 2 symmetry makes easier to solve the equa¬ 
tions because there appears a massless light neutrino. In general, we find all T> tJ parameters 

V?e 2 

depend on the dimensional constant K = Upy- It is true for both normal and inverted 
mass hierarchies. One more interesting result is reached when we take into consideration 
the LFV processes such as Br (/1 —> e + 7) and Br (r —> /i + 7). The current bounds [36] on 
these processes constrain e. For the normal case, one obtains e > 1.97 x 1CT 6 , and for the 
inverted case e > 1.58 x 10 -6 . We also have checked that Y^=i m ^i < 0.23eV, coming from 
Planck [5], and the effective Majorana mass bound m ee < 0.27 — 0.65GeV [37, 38], coming 
from double beta decay experiments, were satisfied. 

After the scalar and neutrino sectors of the model were studied and many of the param¬ 
eters were set, we consider the 7 i R3 (more precisely V DM ) as a DM candidate. We study the 
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bounds coming from the relic density abundance Ddm^ 2 [5] and the direct detection exper¬ 
iments [12-14], Basically, we have worked with three free parameters, k,h 2 , ^ 123 , and the 
DM mass, Mqm- These parameters have been chosen because they play a very important 
role in determining both the IVdm annihilation cross section and the IVdm elastic scattering 
off the nucleon. Roughly speaking, we find that for Mdm 500 GeV, the Ddm ^ 2 < 0.1193 is 
achieved around the resonance regions. For Mdm > 500 GeV there are several regions with 
Odm ^ 2 < 0.1193 aside from resonances regions. This is understood by realizing that the 
couplings of IV D M to scalars (including the Higgs) are proportional to ^p° M , and that for 
Mdm < 200 GeV the main annihilation channels are, in general, mediated by the Higgs. It 
is also observed, that making K 123 bigger, we obtain more regions with Dom ^ 2 < 0.1193 for 
smaller Mdm- This is because K 123 strongly controls the trilincar couplings between scalars 
and thus, the annihilation cross section is larger when K 123 is larger. We have found that 
the relative contributions to the DM annihilation in this model have an intricate pattern. 
It strongly depends on the scalar masses. However, some general conclusions can be drawn. 
For Mdm < 200 GeV the annihilation into gg, W/Z , ll are dominant. For 200 < Mdm < 700 
GeV annihilations into scalars are the most important. Finally, for 700 GeV < Mdm anni¬ 
hilations into N,N t play a important role. 

For DM direct detection, the parameter kh 2 is the most relevant since it is the only one 
which effectively couples IVdm to the quarks in our model. Since nuclei are made of quarks 
(and gluons), this interaction is of supreme importance to the elastic scattering of IVdm off 
nuclei. We found out that if we choose kh 2 = 10~ 2 , our entire curves are below LUX data, 
the most stringent upper bounds on SI DD, however if one chooses a riskier value such as 
kh 2 = 0.1, it can still be lower than LUX, however only above Mdm ~ 500 — 600 GeV or at 
the resonances below that Mdm- 

A final remark concerning to the Z 2 gauge boson is in order. I 11 the region of parameters 
that we have studied, the Z 2 boson does not affect the DM properties. It is because Z 2 is 
heavy since its mass has to satisfy Mz 2 /gB-L > 6TeV [57, 58]. In addition, its mixing angle 
in the neutral current is limited to be tan/3 < 1CU 3 [54-56]. 
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Appendix A: THE MINIMIZATION 

The general minimization conditions coming from dV^-c/dRi = 0, where V^-c is the 
scalar potential in Eq. (2) and R< = {H R} & R , (p 1R , (J) 2 r, 4 >3 r, 4>xr} are the neutral real 


components of the scalar fields, can be written as: 

0 = Vh (2A hVh + k h<s>V§ + KmVfa + ach 2 V(/ 2 + Rh'sV^ + khxV^ x — 2 

— V 2 kh^xV^V^x'j (Al) 

0 = (ac h<s>V r + 2Act.V i/ + + ac$ 2 V(/ 2 + ac$ 3 V^ 3 + k$xV£ x — 2/i|) 

—\[2 k h<i> x Vh V<p x ; (A 2 ) 

o — V<tn ( k hiVh + + 2AiV^ + AC 12 V(/ 2 + AC 13^3 + AClV yf, x — 2/^i) 

+ V < j > 2 V (t>3 (AC123V/) 3 + M23 xV < j > x ) ; (A3) 

0 = ^2 { k H2Vh + + M 2 Vfa + 2A 2 V^ 2 + k 23 ^ 2 3 + « 2 xV^ x — 2/i/) 

+ V^iV03 (Avi 23^>3 + K 123 xV<t>x) I (A 4 ) 

0 = ^>3 { K H3V R + ac$ 3 Vi/ + KwVfa + ac 2 3 V (/ 2 + 2 A 3 VJ 3 + k 3 xV£ x + 3a4xV/) 3 V/ x 

—2/x§) + VfaVfo (2ac 123 V^ 3 + Ki23xV(/>x) > (A5) 

0 = { k hxVh + AC<r.xVi/ + ac ixV^ + ac 2 xV// 2 + AC3xV// 3 + 2\xV£ x — 2fix ) 

—a/2ac h&xV&Vh + Vfo (kusxV^V^ + k' 3 X V^ 3 ) • (A6) 
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